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Abstract. We address the question of how phonons of various frequencies contribute to the superconducting 
transition temperature Tc of a multi-band superconductor by considering the functional derivatives of Tc 
with respect to various intraband and interband electron-phonon coupling functions. A general scheme for 
computing such functional derivatives is developed. The proofs are given that the functional derivatives 
which are diagonal in band indices are linear in phonon energy J7 at small O, while the functional derivatives 
which are off-diagonal in band indices diverge at f2 = as 1/J7. The case of a two-band model for MgB2 
is treated numerically. 

PACS. 74.20.-z Theories and models of superconducting state - 74. 70. Ad Superconducting metals; al- 
loys and binary compounds (including A15, MgB2, etc.) - 74.62.-c Transition temperature variations 
(superconductivity) 
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1 Introduction 

Bergmann and Rainer [1] introduced an important diag- 
nostic tool into the Eliashberg theory of superconductiv- 
ity (for a pedagogical review of Eliashberg theory see [2] ) . 
They considered the functional derivative of the super- 
conducting transition temperature Tc with respect to the 
electron phonon coupling function a'^F{f2) of an isotropic 
(dirty) superconductor 

STr, 



- Tc[a'F{ f2')]) 1^ . 



(1) 



The functional Tc[a'^F{{2)] is defined by the Eliashberg 
equations at Tc [2]. The function STc/Sa^F{n) provides 
answer to the question - How are the phonons of frequency 
i7 effective in contributing to Tc? From the practical point 
of view the functional derivative 5Tc/ Sa^ F{fi) gives the 
change AT^ in transition temperature when a'^F{fi) is 
changed by a small amount Aa^F{fi) (say, by applying 
pressure [3], by alloying [4,5], or by implanting small con- 
centrations of hydrogen into a metal [6]) 



ATr 



STc 



Sa^F{f2) 



Aa^F{n) 



(2) 



The concept of functional derivative 5/5a'^F{Q) was ex- 
tended by Rainer and Bergmann [4] and others to sev- 
eral thermodynamic properties of isotropic superconduc- 
tors (for a review and an extensive list of references see 
[7]). 



The main conclusions of Bergmann and Rainer [1] were 
that 6Tc/ 5a'^F{fl) is always positive (they were able to 
prove this mathematically for the case when the Coulomb 
repulsion parameter ji* [2] is zero), proportional to fi at 
J7 ■^2ttTc (we choose units such that h = 1 and ks = 1) 
and with a maximum at f2 just above 27rTc. Hence, the 
electron coupling to a phonon of any frequency has a posi- 
tive contribution to Tc, but the small values of 5Tc/5a'^F{fI) 
in the low frequency region imply that the changes of 
a^F{n) in this frequency range have no appreciable effect 
on Tc, in contrast to the influence of the low frequency part 
of a^F{fl) on the electron-phonon coupling parameter A 



A = 2 



+ 00 



(3) 
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which is used in McMillan-type interpolation formulae for 
Tc [2]. 

Daams and Carbotte [8] considered the functional deriva- 
tive of Tc with respect to the Fermi surface averaged electron- 
phonon coupling function of an anisotropic superconduc- 
tor and found that 5Tc/5a'^F{Q) diverges at I? = as 
1/fl. Their explicit calculations for a separable model of 
anisotropy a^Fk,k'(^) — (l-l-ak)a^F(/2)(l-l-ak')) with the 
Fermi surface averages (ok) — and (a^) <Cl, showed 
that at small D 6Tc/ da^ F{fl) goes negative and diverges 
as —l/fl. Thus in high purity anisotropic superconduc- 
tors the electron coupling to low frequency phonons de- 
creases Tc, which is analogous to the effect of elastic im- 
purity scattering on transition temperature of anisotropic 
superconductors. However, Daams and Carbotte pointed 
out that for their choice of anisotropy parameter for Pb 
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((a^) = 0.04) 5Tc/5a^F{fl) becomes negative only be- 
low 0.6Tc where a^Fl fl) has very little weight and varies 
as , rendering l/^2-divergence in functional derivative 
harmless. The net effect of anisotropy in the pairing in- 
teraction is to increase the as the virtual scattering 
(k I, — k J.) — > (k' t, — k' i) — > • • • over the Fermi surface 
takes advantage of the regions where the pairing interac- 
tion is large. 

A consensus has emerged (for a review see [9]) that in 
order to describe the superconducting properties of a 40K 
superconductor MgB2 the Eliashberg theory has to be ap- 
plied to a multi-band case with [10] or without [11] gap 
anisotropy on different sheets of the Fermi surface. In [11] 
the Eliashberg equations for an effective two-band model 
of electronic structure and elcctron-phonon coupling in 
MgB2 were solved and the calculated specific heat differ- 
ence between the superconducting and the normal state 
was in good agreement with experiments over a wide tem- 
perature range below Tc- As a model for gap anisotropy 
the two-band model is the opposite extreme to the separa- 
ble anisotropy considered in [8] - there are four Eliashberg 
functions a^FCTCT(/2), a^FT,^{fl!), a'^Fcr^{fi) and a^F^^{Q), 
and, correspondingly, there are four functional derivatives 
of Tc with respect to each one of them (the band off- 
diagonal functions a^F^^^ and o^Ft^^ are related, but are 
different as they are proportional to the partial electronic 
densities of states in tt- and a- bands, respectively). In this 
work we calciilatc the functional derivatives of Tc for the 
two-band model and electron-phonon coupling functions 
presented in [11]. 

The rest of the paper is organized as follows. In Sec- 
tion 2 we present the formalism necessary for computation 
of the functional derivatives of Tc with respect to various 
electron-phonon coupling functions in a multi-band case. 
We also prove that the band-diagonal functional deriva- 
tives 5Tc/6a^Fii{fi), where i is the band index, are pro- 
portional to fi at small i7, while the functional derivatives 
which are off-diagonal in the band indices, STc/Sa^Fij{n) 
with i ^ j, diverge at small fl as 1/J7. In Section 3 we 
present and discuss our numerical results, and in Section 
4 we give a summary. 



2 The functional derivatives 5Tc/5a^Fij{Q) 

In the case of several bands « = 1, 2, . . . with different par- 
tial densities of states Ni the Eliashberg equations for Tc 
do not have the form of a Hermitian eigenvalue problem 
[11]. That is because the interband electron-phonon cou- 
pling functions a'^Fij{Q) and the corresponding Coulomb 
repulsion parameters are proportional to Nj and are 
not symmetric under the exchange of the band indices i 
and J, i 7^ J. As a calculation of the functional deriva- 
tive of Tc relies on the Hellman-Feynman theorem [1], [12] 
which is valid only for Hermitian matrices, it is necessary 
to cast the Eliashberg equations at Tc into a Hermitian 
eigenvalue problem. To this end one first takes the cutoff 
Uc in the Matsubara frequency sums to be large enough 
so that iJ,*j{u)c) = ^ij = V^jNj, where V^j is the Fermi 



surface averaged screened Coulomb matrix element be- 
tween the states in the bands i and j; clearly Vjl = V^'^. 
The electron-phonon coupling functions can be written as 

a^F,j{f2) = a^f,j{Q)Nj with fij{Q) = a^fji{f2). Then 
the Eliashberg equations at Tc take the form 



uJnZi{n) =Wn-\- ttTc ^ Ajj (n - m)Nj 



r+oo 
l)= , 
JO 



Xij{n — m)= I dfia^ fij{Q) 



2f2 



n'^ -h (w„ - Wto)2 



(4) 
(5) 

(6) 



where 4>i (n) and Zi (n) are the pairing self-energy and the 
renormalization function, respectively, at Matsubara fre- 
quency Un = 7rTc(2n — 1) in band i. By defining 



(n) =Mn) / ^\ujrr\Z,{n) (7) 
Eq. (4) takes the form of a Hermitian eigenvalue problem 



i)=eiT)^nT 



\lAn-m)-til 



jm 



\/\0Jn\Zi{n)^\uJm\Zj{m) 



(8) 

where the symmetrized A's and Coulomb repulsion param- 
eters are given by 



XfJn-m) 




^'^ + (W„ - OJmy ' 

(9) 



(10) 



and the eigenvalue e{T) is 1 when T = Tc- In terms of 
A|^(n — m) the renormalization function Zi{n) is given by 
(see Eq. (5)) 



LOnZi{n) =LOn + T^Tc } A* - (n - m) 



jm 




(11) 



Next, in order not to deal with a matrix of unnecessarily 
large size one cuts off the Mutsubara sums in (8) at a 
smaller energy Wg, which is large enough so that Zi(n) ~1 
for |a;„| > Wc, and at the same time rescales /if^- to the 
new cutoff ojc by integrating out the high energy part of 
(t>j (m) as described in [2] . The result is that jil^ in Eq. (8) 
is replaced by iJ.*j{u)c) where the matrix (in band indices) 
fi*{cOc) is related to matrix fi^ by 



A*(wc) 



(12) 



with Ep on the order of the total bandwidth. In the case 
of a two-band model, which we will examine numerically 
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in the next section, the explicit relations between fi*j (uc) 
and iJ^lj are (/i^^ = /x^^) 



s 2x1 Ef 
U)r. 



+ )ln — 



/D,{13) 
(14) 
/D, (15) 



where D is the determinant of i + fi'^ \n{EF/iOc) 
D 



1 + (ML+M:jln^ 



(16) 



Now it is easy to generalize the procedure for calcu- 
lating the functional derivative of Tc given in [12] for a 
single-band isotropic superconductor to the case of sev- 
eral isotropic bands. One finds 



6Tr 



dejTc 
dT 



(17) 



with 



N,. 




2J7 



2J7 



-h (27rT(n -h m - 1))^ 



No 

E 



^f(n) 



/ n-l 

^E 



2n 



/EE<^'fcH- 



The derivative de(rc)/(ir(<0) is conveniently calculated 
in the process of finding the highest T for which the largest 
e(r) in Eq. (8) is equal to 1 and (j)i{n) are the components 
of the corresponding eigenvector of length dNc, where d 
is the number of bands and Nc = [a;c/(27rTc)+0.5], where 
[• • •] denotes the integer part. Note that 0i(n)/-\/w„Zj(n) = 
^/T7i4>i{n) /uJnZi(n) = \fNiAi(n)/ujm where Ai(n) is the 
gap function in band i at Matsubara frequency 

For Q <C 2'kTc one finds from Eq. (18) that the band- 
diagonal functional derivatives are given by 



fc(rc) 

5a^Fu{f2) 



= 2'kTM 



E E 

. i n,m= 

1 



(j)i{n) (j)i{m) 



(27rTc(n - m))^ 



1 V'^nZiin) ^(JmZi{m) 

+ 



(27rTc(n + m - 1))2 
^f(n) 



-E 



^^a;„Z,(n)|J(27rT,0 



AT, 



k n=l 



(19) 



and are linear in i? just like in the one-band isotropic case 
[1]. The band-off-diagonal functional derivatives {i ^ j) in 
the small f2 limit are given by 



Se{Tc) 



27rre 



E 



(f>i{n) 



(j)i{n) 



/EE^^w- 



(20) 



k n=l 



and diverge at /2 = as l/J? just like the functional 
derivative of Tc with respect to the Fermi surface aver- 
aged electron-phonon coupling function of an anisotropic 
superconductor [8]. Note that the sum in the numerator 
of Eq. (20) can be written as A^i X^^Ii -44(71) (zij(n) - 
Ai{n))/LO^ and the sign of 5Tc/ 5a'^ Fij{f2) is determined 
by the relative size of the gaps Ai (n) and Aj (n) near Tc 
in the two bands and by the sign of of Ai (n) near Tc for 
low n since the low-n terms give the largest contribution 
to the sum because of in the denominator. 

The results presented so far apply to a superconductor 
with any number of bands with isotropic intraband and in- 
terband interactions (both electron-phonon and Coulomb). 
In the next section we present numerical results for an ef- 
fective two-band model of MgB2 described in [11]. 



(18) 3 Numerical results for a two-band model 



Following the work of Liu at al [13], Golubov at al [11] 
reduced the four-band electronic structure and electron- 
phonon c;oupling in MgB2 to an effective two-band model 
by exploiting the similarity of the two cylindrical (a-bands) 
and the two three-dimensional (7r-bands) sheets of the 
Fermi surface. We used the a^F^s given in [11] with the 
coupling parameters (see Eq. (3)) X^a = 1.017, A^tt = 
0.446, = 0.212 and = 0.155. The Coulomb repul- 
sion parameters were determined using the ratios of the 
screened Coulomb interaction parameters for MgB2 calcu- 
lated in [14], ^o-a- ■ Mtttt : ^J■a■^^ ■ I^TTo- — 1-75 : 2.04 : 1.61 
: 1.00, the density of states ratio N^/Nc, consistent with 
the ratio AcrTr/A^-o- for the spectra that we used, and us- 
ing equations (10) and (13-16) with Ep set equal to the 
TT-bandwidth of 15 eV [15] and uJc — 0.5 eV. These con- 
straints leave the single fitting parameter Hcra which was 
fitted to the experimental transition temperature of 39.4 
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Fig. 1. The band-diagonal functional derivatives and rescaled ^ ^^^^ ^-^ ^ ^^^^ ^ ^^^^^^^ ^^^^^^ ^^^^^ 

electron-phonon coupling functions over the entire phonon en- ^^^^ ^^^^^ oia^F^^iQ) and a^F^^{n) that are near the 



ergy range in MgB2 



maxima in the corresponding functional derivatives are drawn 
to scale. 



K. The results of the fit were: Ai*<,(wc) = 0.19627, ti%^{uJc) 
= 0.19561 and n%^{ujc) = M^^K) = 0.04948. 

The calculated functional derivatives are shown in Fig- 
ures 1-3. In Fig. 1 we present the band-diagonal functional 
derivatives 5Tc/5Q^Faa{^)i 5Tc/5a'^FT^.,r{fi) and rescaled 
o?Fcrcr{^) and a^i^7r7r(^) ovci the entire phonon energy 
range in MgB2. The spectra a^i^(jo-(i7) and a^F.^T^{fl) are 
scaled down by a factor of 20 so that their shapes over the 
entire range of phonon energies are seen on the scale set by 
the size of the corresponding functional derivatives. This 
leaves parts of a^F^^l fi) and Oi^ F-j^T^i^f^) that are near the 
maxima in corresponding functional derivatives invisible, 
and in Fig. 2 we redraw a part of Fig. 1 over a smaller 
energy range with the electron-phonon coupling functions 
shown to scale. As we proved in Sect.|21 5Tc/ Sa^ F„„{fl) 
and 5Tc/ Sa^ FT^T^{fl) are linear in Q at small energies and 
vanish at = 0. Moreover, both functional derivatives are 
positive with broad maxima at 9-10 times fesTc, similar 
to what one finds in the isotropic single-band case [1],[12]. 
The difference in sizes of the two functional derivatives in 
Figs. 1 and 2 is due to the difference in sizes of the gap- 
functions Aa{n) and AT^{n) near Tc in the two bands as 
can be deduced from Eq. (18) (or Eq. (19)) by noting that 
the normalization factor X]n('^o-("') + '^^("')) the same 
for both functional derivatives and it's size is largely deter- 
mined by 4>ain), which is larger than (t>-K(n) (^cr(l)/07r(l) 
= 3.4). Note that in the isotropic one-band case the scale 
of the functional derivative of Tc varies roughly as l/(l-f A) 



0.6 



0.2 







-0.2 





|i 
|i 

J^) ^\\ 


- 1 a¥ 


lili 




ii il 
il ) 


^^^^ 


1^ Ii 11 
\\f 'ii M 


; 

// 
// 


1 litt — ;^-^,-^!^_ ■ 
l\ /; V '■■■■■■...\J 
\</ ■■■NJl 









50 100 

J^(meV) 

Fig. 3. The band-off-diagonal functional derivatives and 
electron-phonon coupling functions for MgB2. 
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[16] and from Fig. 1 it is clear that such a "rule" cannot 
be applied in determining the relative size of the band- 
diagonal functional derivatives of Tc in a multi-band case. 

In Fig. 3 we present the band-off-diagonal functional 
derivatives JTc/fo^i^crTrl^) and 5Tc/(5a^F^o-(^) together 
with the corresponding electron-phonon coupling functions 
As we proved in Sect. [21 both of these functional deriva- 
tives diverge at = as but they also have opposite 
signs. Both a'^Fc-,r{fl) and a^F-^cri^Q) vary as fl^ in the 
limit — > so that l/X?-divergences in the corresponding 
functional derivatives are integrable, as can be deduced 
from Eq. (2). 

The difference in signs between 6Tc/6a^FcrTr{f2) and 
STc/So^Ft^^^Q) is related to the fact that near Tc A^{n) > 
AT^{n) with both gap functions positive at low n (see 
Eq. (20) and the subsequent discussion in Sect.|21). In order 
to illustrate the importance of the sign of the smaller gap 
At^^h) at low n close to Tc, we have computed the func- 
tional derivatives for the case when a^FT^Tri^), oi'^F<j-k{^) 
and c?FT^i,(fl) were scaled down by a factor of 10, with 
o?-F„a{n) and M^a(wc), Ai^^(wc), i^lA^c), ^^lai^c) left un- 
changed. This produced negative Atttt — fJ-'^Tr, ^-rra — ^J'*^^a 
and Ao-TT — Mcttt which resulted in a solution where A„ {n) 
and Zi^(n) near Tc have opposite signs at low n. The 
corresponding functional derivatives are shown in Figure 
4. Now, both band-off-diagonal functional derivatives are 
negative. Note that the scale of STc/So^Ft^.^ is roughly two 
orders of magnitude smaller than the other three func- 
tional derivatives. 

It is important to point out that the calculated Tc 
for the parameters in Fig. 4 was 43.7 K - substantially 
higher (by 10%) than the Tc of 39.4 K obtained for the 
parameters used to produce the results given in Figs. 1-3. 
In fact, even higher Tc of 45.1 K was obtained by set- 
ting a^F^^, a^FcTTT and a^F^^ identically equal to with 
a^Fccr and all the Coulomb repulsion parameters left the 
same - i.e. no attractive interaction in tt-tt and ct-tt chan- 
nels! The corresponding functional derivatives were sim- 
ilar to those shown in Fig. 4 and they explain why Tc 
is reduced as the couplings a^F^^, a^FcT^ and o^Fct^ are 
turned on from 0: bTcj ^o?FaT^ and bTd ^o^Ft^ci are both 
negative and much bigger in absolute value than than the 
positive STc/So^Ft^t^. We finally note that if all interac- 
tions, both electron-phonon and Coulomb, in tt-tt and ct-tt 
channels are set equal to zero, which effectively reduces the 
the two-band model to one-band model, the calculated Tc 
was 44.6 K. This is a half degree lower than what was 
obtained by turning off only electron-phonon interaction 
in TT-TT and ct-tt channels. 



4 Summary 

We have developed the general formalism for calculating 
the functional derivatives of Tc with respect to electron- 
phonon coupling functions for a superconductor with sev- 
eral bands with isotropic intraband and interband inter- 
actions (electron-phonon and Coulomb) . We proved rigor- 
ously that the band-diagonal functional derivatives 5Tc/ 
Sa^Fii(f2) are linear in fi at small f2, as in the isotropic 
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Fig. 4. The functional derivatives for the case when the 
strengths of a^F^^, a^Fcrrr and a^F^a- were scaled down by 
a factor of 10 with a^Fa-a- and /i*((^c)'s left unchanged com- 
pared to those used in Figs. 1-3. 

single band case [1] . At the same time we proved that the 
functional derivatives which are off-diagonal in the band 
indices, STc/5a^Fij{f2) with i ^ j, diverge at small SI as 
The calculation was carried out for a two-band model 
of MgB2 using the electron-phonon coupling spectra given 
in [11] and the ratios of the screened Coulomb interaction 
parameters given in [14]. We found that the band-diagonal 
functional derivatives are both positive with broad max- 
ima in the range of 9-10 times fc^Tc, similar to the sin- 
gle band isotropic case. However, the functional deriva- 
tive with respect to intraband electron-phonon coupling 
function for the band with the smaller gap (tt band) was 
found to be much smaller than the corresponding func- 
tional derivative for the band with the larger gap (cr band) . 
The functional derivatives with respect to the interband 
electron-phonon coupling functions were found to diverge 
as l/!7 at 0, but had opposite signs over the entire 
range of phonon energies for the parameters given in [11]. 
We found that the signs of these off-diagonal functional 
derivatives are determined by the relative signs of the gap 
functions near Tc at low Matsubara frequencies in the two 
bands and that, in general, it is possible to have both 
band-off-diagonal functional derivatives negative. 

The results found here give a better insight into the 
questions - What is the effect of phonons of frequency S2 
on Tc through their couplings to electrons via various band 
channels? Are all coupling always contributing positively 
to Tc, or are some of them pair-breaking? The answers to 
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these questions are provided in Figures 1-3. 
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